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The well-known phase structure of the two-dimensional sine-Gordon model is recon- 
^ . structed by means of its renormalization group flow, the study of the sensitivity of the 

dynamics on microscopic parameters. Such an analysis resolves the apparent contradiction 
between the phase structure and the triviality of the effective potential in either phases, pro- 
vides a case where usual classification of operators based on the linearization of the scaling 
relation around a fixed point is not available and shows that the Maxwell-cut generates an 
unusually strong universality at long distances. Possible analogies with four-dimensional 
Yang-Mills theories are mentioned, too. 
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Introduction The phase structure and renormalizability of field theoretical models 
are well understood in the context of the renormalization group: Phase transitions 
belong to the singularities of the effective IR coupling constants as the functions of 
the UV parameters and the renormalizable parameters are the relevant or marginal 
coupling constants of the UV fixed point. The two-dimensional sine-Gordon (SG) 
model has a well-known phase structure and renormalization group flow but does 
not easily fit into the general scheme. Our aim in this work is the clarification of 
these issues by a careful renormalization group study of the SG model. 



The SG model, defined by the action 
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has an ionized (massless, strong coupling, non-renormalizable) phase for P"^ > Stt 
and a molecular (massive, weak coupling, renormalizable) phase for < 8tt. The 
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molecular phase^s^erturbatively equivalent with the neutral sector of the massive 
Thirring model [T, 2] and the neutral Coulomb-gas |5] • Perturbation expansion was 
used to conjecture the same universal behavior around = Svr ||4i [SD and phase 
structure [6] as that of the planar X-Y model. Simple comparison of the lattice 
regulated SG model with the planar X-Y model provides a non-perturbative renor- 
malization group flow for the SG model [7]. The direct lattice regularization of the 
SG model makes it evident that our model, given by Eq. ©, has aperiodic kinetic 
energy. In fact, the first term in the action is periodic for space-time independent, 
p = Fourier mode only. Its aperiodicity for space-time dependent fluctuations 
suppresses configurations which correspond to the vortices of the X-Y model. As 
a result, the model has no or few renormalizable coupling constants in the ionized 
and the molecular phases, respectively [170. It has been argued that the appropriate 
order parameter, provided by the soliton dynamics, is the topological susceptibility 
because the symmetry with respect to the fundamental group 7ri(f/(l)), 

<Px ^ <Px + A, (2) 

where 

27r 

A = j (3) 

is the period length of the local potential, is broken dynamically in the molecular 
phase. 

But some complication arises by the similarity of the deep infrared IR dynamics in 
the two phases of the SG model [8]. The rather trivial observation, namely that the 
only non-singular effective potential which is periodic and convex is the constant 
shows that the Maxwell-construction washes away the differences of the phases in 
the deep IR regime. How can this be reconciled with the overwhelming evidences 
about the two differing phases and the different numbers of the renormalizable 
parameters? We show that the difference of the two phases appears in the effective 
potential expressed in units of the running cutoff. This potential produces weak 
effects in the deep IR but can be used to identify the phase structure because it 
produces singularities in the sensitivity of the dynamics at a given observational 
scale as the function of the bare parameters, given at the cutoff scale [9]. Such a 
global use of the renormalized trajectory is needed because of the non-triviality 
of the IR scaling laws, the inherent non-linear nature of the renormalization group 
trajectory at the IR fixed point of the ionized phase. 

Blocking and condensate The blocking relation for the Euclidean Wilsonian action, 
Sk [0] to be integrated from the initial condition imposed at the cutoff /c = A, is 

g-Sfe_AfcM ^ f /}U']e-5fe[<^+'^'l (4) 



where the Fourier amplitudes of and (f)' are restricted to wave numbers < < 
k — Ak, and k — Ak < \p\ < k^ respectively [10], in the spirit of the differential 



renormalization group (RG) |lllj,[l20. The functional integral is over modes whose 
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scale belongs to a finite interval, therefore the dependence of the blocked action 
Sk~Ak[<P] on the parameters of the original action, S'jt[0], is regular flB*], no phase 
transition-like singularities flT4l ITsl 16 1 are allowed. The loop-expansion of the 
functional integral gives 



Sk-Ak[4>] = Sk[(f> + 0o] + -Tr'ln 



K.nJ"Sk[<l> + <P'o] 



(5) 



where 0g is the background field dependent saddle point nln\ 



The local potential approximation (LPA) for the SG model consists of the projec- 
tion of the evolution equation into the functional space 



Sk 



l{dM' + Uk( 



(6) 



where the local potential is given by the Fourier series 

oo 



(7) 



n=l 



The Wegner-Houghton equation [lOd reads 

(2 + kdk) Uki(t>) = In (l + f/f (0)) , 



(8) 



when the vanishing of the saddle point is assumed, Uk = k "^Uk is the dimension- 
less potential and ui^\(l)) = d'^Uki'p). 

The saddle point which indicates the presence of a condensate in the system, the 
hallmark of spontaneous symmetry breaking, appears at the scale where the Eu- 



clidean inverse propagator, 6 Sk 



, looses its positive definiteness. The neg- 



ative value of the inverse propagator indicates that the action has unstable mini- 
mum for vanishing strength of the field components to be eliminated. The action 
is bounded from below therefore the true minimum is degenerate, implying highly 
populated modes, a condensate whose strength is stabilized by the balance between 
kinetic and potential energies. 

The loss of positive definiteness of the propagator corresponds to the change of 
sign of the argument of the logarithmic function in Eq. ([8]) in our approximation. 
Similar phenomenon has already been observed in the case of the scalar 0"^ model 
in Ref. [18], too. Once the tree-level contribution to the blocking relation appeared 
we neglect the loop contributions and the evolution is governed by the tree-level 
blocking relation Sk-Ak[<P] = niin(^/(S'fc[0 + 0']). We search the saddle point among 
plane waves, 0^ = pk cos{knk ■ x + 6k) with Uk being a unit vector, and the result 
isfli 

mjn + 2 y duUk{(j) + 2pcos{7iu)) , (9) 



Uk~Ak{4>) 
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replacing Eq. ([8]) in the unstable region. It is easy to show that the RG flow in LPA 
preserves the period length of the potential. 

The saddle point reflects the macroscopic population of a mode and is a well known 
phenomenon in theories with spontaneously broken symmetry. For instance, the 
symmetry broken vacuum of the simple (f)^ model has a condensate in the mode 
p = 0. What is new in the saddle point occurring at finite k in the blocking equa- 
tion is that it is space-time dependent. There are such inhomogeneous, macroscop- 
ically populated modes in the (f)'^ model if the absolute magnitude of the space-time 
average of the field expectation value is constrained to smaller values than in the 
symmetry broken vacuum. Such a system is in a mixed phase because the constraint 
is satisfied by a domain structure where each domain has a stable, clusterising vac- 
uum expectation value of the field. The inhomogeneous saddle points represent the 
domain walls. 



An important consequence of the inhomogeneous saddle point is the appearance of 
soft modes, the analogies of the Goldstone modes, even if the internal symmetry, 
broken spontaneously, is discrete. In fact, the parameters 9^ and are zero modes 
arising from the breaking of the continuous space-time symmetries, translations 
and rotation, respectively, by the inhomogeneous saddle point. They are the soft, 
non-perturbative modes of the domain structure and generate the Maxwell-cut, the 
degeneracy of the true ground state as the function of the field expectation value 
within the mixed phase HItIi . 

Do we expect condensate in the SO model? The stability of the kinks originate 
form the spontaneous breakdown of the fundamental group symmetry Thus 
condensate is expected when so long the solitons are stable. The presence of virtual 
kink-anti kink polarisation of the vacuum indicates that this condensate contains 
inhomogeneous modes, too. 



After having clarified the dynamical role played by inhomogeneous condensates in 
our model we return to the evolution equation ([5]). The change of the sign of the 
propagator signals a singularity in the second term on the right hand side, an IR 
Landau pole which is always an artifact of the truncation. Although the conden- 
sate modifies the dynamics and restores the positive definiteness of the propagator 
nevertheless the question remains: What happens with the exact solution to the 
evolution equation? An important test is passed by the numerical integration of the 
evolution equation, it reproduces regular scale dependence in the blocking relations, 
as expected. But the sharp cutoff in the evolution equation (U) for the Wilsonian ac- 
tion prevents us to go beyond the LR^. One should use smooth cutoff which can be 
implemented for the Wilsonian action [11] or the effective action il2ll . The former 
would be more appropriate in our case because the evolution equation for the effec- 
tive action leads to degenerate action [19] but this is the natural manifestation of the 
convexity of the Legendre transform and completely hides the dynamics behind the 
Maxwell-cut. We do not pursue this direction of improvement in the present work. 
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instead we assume that the appearance of inhomogeneous saddle points in our trun- 
cation scheme is already the correct reflection of macroscopically populated states, 
a spontaneously broken symmetry with respect to the fundamental group transfor- 
mation in the vacuum. 

The solution of the RG equation ([8]) in the UV scaling regime, fc^ 3> | uf'' \ , is 



after ignoring contributions 0{\Uj^ \'^/k'^) in Eq. ([8]), displaying the well-known 
critical point at = Stt. The relevant coupling constants of the UV fixed point 
are called renormalizable because those parameterize the dynamics at finite scales 
when the cutoff is removed. The two-dimensional scalar model with polynomial 
interactions possesses infinitely many renormalizable operators. In contrast, the pe- 
riodicity of the potential of the SG model seems to make all parameters of the local 
potential UV irrelevant, i.e. non-renormalizable in the ionized phase, > Svr and 
to allow only a few renormalizable parameters in the molecular phase, < Svr. 
But one has to consider the renormalized trajectory globally, by taking into account 
the IR scaling laws in order to find the free parameters of the renormalized dynam- 
ics. In fact, the number of these parameters might be more or less if the IR scaling 
has new relevant operators or stronger universal features, respectively [@]. We now 
present a global study of the renormalization group flow to point out that the local 
analysis at a given fixed point is actually not reliable in the SG model. 

The solution of the RG equation can only be obtained numerically in the IR region 
because of the mixingof the different Fourier modes and the appearance of con- 
densate for (3'^ < Svr [jsl]. It was also shown in Ref. [8] that the effective potential 
of the SG model is rendered to a field-independent constant by the requirements 
of periodicity in the internal space and convexity. Therefore, all the dimensionful 
couplings Unik) should tend to zero in the IR limit. Nevertheless according to the 
expectation value of the topological susceptibility as the disorder parameter of the 
model one can distinguish two phases, an ionized one for > Svr, and a molecular 
one for < Stt. As we show below this feature is reflected in the different shape 
of the dimensionless effective potential in the two phases. 

Ionized phase There is no spinodal instability and the renormalized trajectory can 
be well approximated in the IR scaling regime by the simple power law 



with r] > 0,kQ being some scale parameter. In fact, the evolution equation reads as 




(10) 



Un = Cn{k/ko) 



(11) 



(2 + nr])ncnk' 




■n'^Cnk^'^ + -13^ ^^n,s{2 + 57/)^^" (12) 



1 



s=l 
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with An,sik) = in 
n = 1 one finds i] 



to the recursion relation 



— {n + s^Un+s, when this assumption is made. For 

2 > 0, Ci = 'Ui(A)(A;o/A)'' and the cases n > 1 lead 



n{2 + nri -ri^^) 



(13) 



expressing c„ in terms of ci, c„ = (— 1)"+^m"(A)_R„, where Ri = 1 and the Rn 
satisfying the recursion relation (fT3l) becomes independent of the bare couplings. 
This scaling law is confirmed by the numerical results, shown in Fig. [U for r/ = 1, 
i.e. (3"^ = VItx. Further numerical support is that the ratio i?„ = \un{k)\ / u^{k) = 
|cn|/(ci)" becomes A;-independent and the local potential becomes independent of 
M„(A), 7T, > 1 in the IR region. 




k/A 

Fig. 1. The IR scaling law at = 12tt for various Fourier amplitudes. 

A remarkable complication takes place at the IR fixed point, tt* = 0. The upper 
harmonics with n > 1 decrease too fast as — > 0, 0{un) = 0{ui), and the 
renormalizable trajectory is not linearizable around the fixed point. The problem is 
that one tacitly assumes in the usual argument in linearizing the blocking relation 



u' - ut 



- u'-JdmBniu") + 0{{U'„ 



(14) 



that the deviation of each coupling constant from the fixed point values is of the 
same order of magnitude. As a result, we do not have the usual classification of 
operators and a more complete study, based on the global features of the renormal- 
ization group trajectory is needed to determine the number of free parameters of the 
renormalized dynamics. The results like the one shown in Fig. [T] indicate that scal- 
ing laws characterized by critical exponents are actually recovered in a non-trivial, 
non-linear manner with no relevant parameter. This result is in agreement with the 
picture suggested by the analogy with the X-Y model [7] where the vortex fugacity 
is found to be the only relevant parameter in this phase. But these configurations 
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have divergent action for the SG model defined in terms of non-periodic space- 
time derivatives. In other words, our results correspond to the vanishing fugacity 
hyperplane in the space of coupling constants. 

Molecular phase The renormalized trajectory follows scaling laws similar to those 
of the ionized phase at the beginning. Namely, the asymptotic scaling law (flOl ) in 
the UV regime ends at a crossover beyond which the scale dependence m„ ~ /c"'' is 
encountered. The crucial difference between the two phases appears in this region 
in what this scaling law does not extend to A; = as in the ionized phase, rather 
it is interrupted by a further crossover. The period length of the potential is larger 
than in the ionized phase and as a result the symmetry 4>x ^ 'Px + Svr//?, belong- 
ing to the fundamental group of the theory, is broken spontaneously at low energy 
where the potential has more chance to localize the field around a given minimum 
by the formation of inhomogeneous saddle-points to the blocking relation. The sec- 
ond crossover mentioned above, a spinodal instability appears when the propagator 
diverges, fcgj + f/fcgj(0) = 0. The blocked action can easily be determined when 
(3'^ — > Stt from below because Si is the only renormalizable coupling constant. Due 
to fcsi/A ^ 1 the UV irrelevant couplings die out and it is sufficient to keep track 
of ui only, i.e. one can estimate the value of fcsi by calculating the potential Uk{(t>) 
as if it would contain a single Fourier mode only. Using Ui{k) = Ui{A){k/ A)^ we 
find 

kl, = A'{(3'MA))^ . (15) 

We used here the fact that the minimum of U'^{(j)) lies at = in the case of a single 
coupling. The numerically determined running of the coupling constant ui{k) is 
depicted in Fig. [21 Note that the length scale l/A^si where the spinodal instability 
occurs is not an analytic function of /? at = Svr. This situation resembles to the 
non-analyticity of the correlation length at the Kosterlitz-Thouless phase transition 
point 




10'^ lO" 10"'' 10"^ 10'^ lo' 



k/A 

Fig. 2. Scale-dependence of the couplings tt„, with n = 1 . . . 4 at = 47r. The scale A;si 
where the spinodal instability appears is shown explicitly. 
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Once the cutoff has reached the spinodal instability region its further decrease 
strengthens the condensate and flattens the potential in order to arrive at a con- 
stant potential for k ^ 0, being the only function which is simultaneously convex 
and periodic [SJ. The numerical results, displaying this phenomenon are presented 
in Fig.[3l 




Fig. 3. The dimensionful potential and the amplitude of the condensate are plotted for de- 
creasing values of the cut-off for = An as the functions of the homogeneous background 
field. The potential shown by solid line corresponds tok = k^i ~ 0.0046 when the spinodal 
instability appears. The potential flattens out and the condensate grows as k is decreased. 



Though the dimensionful effective potential is flat in both phases, the dimensionless 
effective potential f/o(0) is not flat in the molecular phase where one has Uo{(p) = 
-|((/) - nA)^ for {n - i)A <</><(« + |)A, with n = 0, 1, ... and A being 
given by Eq. (|3]). Both the linear dependence of the strength of the saddle point on 
the background field and the quadratic shape of the dimensionless potential can be 
understood by simple analytical considerations bill . 



Note that the spinodal instability and condensation also induce an effective poten- 
tial of parabolic shape when the periodic bare potential is replaced by a polynomial 
one. The only difference is the lack of periodicity in the latter case. The universal 
parabolic shape is a direct consequence of the LPA and the Maxwell-cut which 
have to be performed on the unstable region of the effective potential lll7l l20ll. 



In order to make the approach of the potential to a parabola more apparent we 
numerically determined /9^m„(0) and found that it is a universal function of n, 
P^UniO) = 2(-l)"/n2 for mi(A) > 0, cf the inset of Fig. SI The corresponding 
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Fig. 4. The magnitudes of the first few dimensionless couplings are plotted as the function 
of /3. The inset shows that the /3-dependence can be easily factorized and the index-depen- 
dence also obtained. 



dimensionless effective potential is 

P ^1 ^ 



(16) 



apart from a field independent constant. This parabola is repeated periodically 
along the (j) axis. It is apparent that such a periodic potential is the solution of 
Eq.®. 



How is the parabolic shape formed as we approach the IR end point? The exactly 
parabolic shape part of the dimensionless blocked potential appears at = k^i 
only and spreads over larger field region as the cutoff is further lowered. But the 
potential approaches a parabola- looking shape already before the appearance of the 
condensate as demonstrated in Fig. [51 where the difference Af/fc(0„) = Uki<Pn) — 
(— is plotted at 0„ = nvr/ (10/5) with n = 0, 1, . . . , 10. The gradual approach 
of the potential to a parabolic shape is a precursor of the condensation. According 
to the numerical results Af/fc(0„) follows a power law behavior in a rather small 
region above ksi. The exponent characterizing this power law depends on P, the 
value u = 3.75 was found at P'^ = Atx with u decreasing as [3 is increased. Once the 
parabolic shape is approximately installed for k slightly above fcsi the potential does 
not suffer a sudden change with the appearance of the condensate and it assumes 
its parabolic shape very rapidly below /csi in the whole period length. 

The insertion of an aperiodic parabola section into a periodic potential and the join- 
ing of functions with different analytical structures are rather involved issues. The 
deviation of the dimensionless potential from the parabolic form around the match- 
ing points are found to be weak but not sufficiently fast converging in the Fourier 
expansion. The increase of the order of the truncation of the Fourier series makes 
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k/A 

Fig. 5. The deviation of the dimensionless blocked potential from the parabolic shape at 
different values, (pn = mr/{W(3) with n = 0, 1, . . . , 10, of the field variable. The deviation 
increases with increasing n. 



the spread of the parabolic shape faster at k ^ kis but the second derivative of the 
potential increases at the matching points without bound, too. Unfortunately the 
Wegner-Houghton method does not allow us to go beyond the LPA where devia- 
tions from the parabolic potential are compatible with the Maxwell-cut. 

The dimensionless potential (fT6l) displays a non-trivial structure in the molecular 
phase and is independent of the bare coupling constants, either relevant or irrelevant 
at the UV fixed point. Such a surprisingly strong super-universality renders the deep 
IR physics completely parameter-free. Though the UV relevant couplings influence 
the dynamics at finite scales, the physics becomes parameter-free well below the 
condensate scale. 



The dimensionless local potential reveals that there are different dynamical mecha- 
nisms leading to the Maxwell-cut and distinguishes the two phases of the SG model 
but this difference is vanishing in fixed, dimensionful units. Nevertheless the phases 
can still be distinguished by means of the renormalization group flow as shown in 
Fig. [6l It remains to be seen if this difference rests qualitatively valid when one 
goes beyond the LPA. In particular, it would be interesting to decide whether the 
couplings Un{k 0) remain finite |22] or diverge flj] beyond the LPA when the 
parameter /? exhibits also scale-dependence. 



Sensitivity matrix The sensitivity matrix is the systematical tool to distinguish 
phases of a theory according to the global features of the renormalization group 
flow. Let us consider the bare coupling constants, Un{k) = -E„(tt(A), r), as a func- 
tion of the ratio of the observational scale and the cutoff, r = k/ A and the bare 
parameters, identified by the initial condition m„(A) = En{u, 1). The sensitivity 
of the renormalization group flow at scale k to an infinitesimal change of the bare 
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Fig. 6. The phase structure of the sine-Gordon model. The dotted line corresponds to 
tti(O) = const. //3^. The critical value = Svr separates the two phases of the SG model. 

couplings is characterized by the sensitivity matrix |@] 

. ?M^. (17) 

Phase transitions can be identified with the singularities in Sn,mir) as the function 
of the parameters of the theory in the limit r 0. The singularities developing as 
r ^ with fixed k and A ^ oo correspond to quantum phase transitions because 
they are driven by short distance phenomena. The traditional phase transitions, in- 
duced by the low energy, long range modes belong to singularities realized with 
fixed A and A; ^ 0. 



The UV scaling laws are given by Eq. (flOl) . 

5„,™(r) = 5„,„r^("'^'-8-) (18) 

in either phase. The situation changes significantly in the IR region, r <^ 1, where 
the scaling law 

u^{k) = {-l)"+^RnuUk) = (-l)"+^i?„<(A)r"i^(^'-8") (19) 

of the ionized phase yields 

Sn,mir) = 5™,i(-l)"+'i?„nMri(A)r"^(^'-«-). (20) 

The same IR scaling behavior is realized in the molecular phase for k > ksi as 
in Eq. (fTTI) . c.f. Fig. [21 thus the Eq. (|20|) remains valid when < Svr. The low 
energy dynamics looses sensitivity on any bare coupling constant when i] > 0. 
This is supported by the numerical results for the sensitivity matrix, shown in Fig. 
|7] with dashed line. The lack of sensitivity is established in two different ways: 
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both in the UV (A — > oo) and the IR (k — > 0) limits. The renormalized SG model 
becomes trivial and has no free parameter in this phase even if the problem with 
the linearization of the blocking relation around the IR fixed point prevents us to 
classify the operators in the usual fashion. 




r 



Fig. 7. The numerically calculated sensitivity matrix at Un(A) = 10 ^5n,i- The slope in 
the ionized phase is in an excellent agreement with t] = 1. 

The sensitivity matrix, given by Eq. (l20l) . has a singularity as the function of the 
coupling constants at = Stt when the limit r — > is made due to scaling laws for 
k > ksi, in agreement with the UV origin of the phase transition in the SG model, 
c.f. the discussion around Eq. (fT5l) . The renormalized dynamics at such a scale k 
has indeed few free parameters, the relevant coupling constants of the UV scaling 
law. This phase transition is shown by the continuous line in Fig. |7]for r > 10^^. 
But the approach of the potential to a parabolic shape as ~ fcsi stops this trend. 
The dynamics depends less and less on the high energy, short-distance parameters 
as k is decreased below ksi and the super-universality developing in this regime 
renders the dynamics free of any parameter as A; ^ 0. 

The dimensionful potential shows the same, flat shape in both phases. Despite the 
clear difference of the dimensionless potential and the sensitivity matrix in the deep 
IR regime, indicating that the dynamics, observed in units of the low cutoff differs 
significantly in the two phases, the IR dynamics of these phases do not look dif- 
ferent when expressed in fixed units. It is easy to understand this apparent contra- 
diction. The two phases are similar in their IR plane-wave dynamics, both decou- 
ple them. But it is well-known that the SG model has other important degrees of 
freedom, solitons, due to the periodicity of the potential. Their finite size or mass 
introduces another important scale where the dynamics differ in the two phases. 

Summary The phase structure and the number of free parameters of the SG model 
were studied in this work. The phase structure was identified without using topo- 
logical disorder parameters, by means of the renormalization group flow. We had to 
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rely on the global analysis of the renormalization group flow due to an unexpected 
problem with the linearization around the IR fixed point of the ionized phase. Fur- 
thermore, it has been clarified that the sensitivity of the dynamics on the renormal- 
izable coupling constants is washed away and a super-universality is generated in 
the IR limit of the molecular phase by the Maxwell-cut. 

These features of the SG model show possible interesting analogies with four- 
dimensional Yang-Mills theories. Both models have periodic variables. The peri- 
odicity of the gauge model, the fundamental group symmetry is supposed to play 
a key role in establishing the confining forces between charges. The renormaliz- 
able, asymptotically free coupling constant of the Yang-Mills theory and the SG 
model in the molecular phase increase with the observational length scale. Such 
an increase generates condensate and spinodal instabilities in the vacuum of both 
theories. The speed of sound is vanishing in the mixed phase with spinodal insta- 
bility, indicating the emergence of color confinement, the absence of asymptotical 



charged plane- wave states [|20n . The exact scattering matrix constructed in terms 
of soliton-anti soliton bound states of the SG model suggests that the asymptotic 
states of this two dimensional model are made up by these composite bound states. 
The SG model realizes a mechanism to remove any free parameter in the dynamics 
of the condensate, well below its energy scale, in a manner similar to the long range 
dynamics of the Yang-Mills vacuum. 
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